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METHOD  OF  APPROXIMATE  CALCULATION  OF  AERODYNAMIC  CHARACTERISTICS  OF 
CAMBERED  WINGS  WITH  VARIABLE  SWEEP 

A.  I.  Pastukhov,  Department  of  Hydromechanics 

The  wing  is  replaced  by  a system  of  attached  vortices  of 
variable  intensity  ever  the  span,  distributed  continuously  over  its 
middle  surface.  Emerging  from  each  point  on  the  vortex  surface  are 
free  rectilinear  vortices  with  variable  angles  cf  decalage  over  the 
chords  and  span  and  arranged  in  the  vertical  plane  parallel  to  the 
root  section. 

In  calculating  inductive  velocities  the  system  of  attached 
vortices  is  assumed  to  be  equivalent  to  the  vortex  system  of  the 
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projection  of  the  Biddle  surface  onto  the  plane  which  passes  through 
its  root  chord  in  a direction  parallel  to  the  span. 

The  angle  of  sweep  of  the  i-th  portion  of  the  attached  vortex 
(Pig.  1)  is  deterained  by  the  expression 

W2)  = **  C . 

where 

f Cm./  C • fp  ? v <»  *,.)  • 

1=-^-  is  ths  relative  coordinate  of  an  elementary  lifting  vortex  in 
the  root  section. 


Assuming  that  the  vortex  density  f(zit>°0  at  section  point  and 
in  the  root  section  J (*,«<)  are  related  by  distributing  function 
in  the  fora  of 


(i) 


for  the  intensity  of  the  elementary  lifting  and  free  vortices  we  can 
get,  accordingly. 


ft*.*)  ?(*„.«)  ir  «»*,** 
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“here  local  narrowing  of  the  wing  projection  is  determined  by 
expression  s 


For  calculating  inductive  velocities  the  v-th  seguent  of  the 
apes  is  broken  down  into  A,  bands  of  width  L*/n>.  Quantities 

assuaed  to  be  constant  and  equal  to  their  values 
in  the  aiddle  if  of  the  bands.  In  this  case,  after  switching  to 
angular  coordinate 


i = -coj8  ; i a . coj8 


for  the  nornal  velocity  conponent  of  the  downwash,  caused  by  the  free 
vortices  at  the  point  with  coordinates 


j'  =-acoj8'; 
<■ 

we  get,  approxinat ely : 


t- 


* /L** « - 4 j Tjlt - '(f£5^rr'Hh’ 


h 

it 


(2) 


v 1 p*  <9  *:  - <9  <>(&,)  * * -fr.  *;*(£■') * ; • 


(3) 

<•) 


where 


1 
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The  velocity  component  of  the  downwash  W**  is  calculated  by 
(2)- (4),  replacing  cm ^ in  (2)  by  J»n/^  ; the  component  is 

calculated  by  replacing  by , uhere 


r (y  )**  #,ji"Vr,(ry (l* 


r(r)s  IV IV  (vv)^*(vv)^ 

By  using  the  sane  assumptions  for  the  velocity  caused  by  the 
attached  vortices,  ve  can  get 


• here 


Vl*,*r'(T')']r 

In  conpiling  the  integral  equation  the  velocities  which  are 
normal  to  the  sing  projection  are  considered  to  be  directed  along  the 
local  noraals  ( nf<  in  Pig.  2)  to  the  niddle  surface.  Then,  the 
condition  of  iapsrnenbility  is  written  in  the  fern  of 

_5««_  f Tc,  _ y cff,y  - ZiL  fir  jinsdB , »5» 

/*  J coj»-cm*'  • ’*  r«  **  J1"  * 151 


where 


<.(•>)*« 
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1*  angle  of  attack  of  the  plane  of  the  wing 

projection; 

” angle  of  twist  of  the  Diddle  section  of  the 

% ~th  plane  in  relation  to  the  root  section; 

(•')  - angle  of  curvature  of  middle  line  of  section; 

- angle  between  local  noraal  and  plane 


If  we  seek  the  solution  to  (5)  in  the  form  of  a series  [3],  [4], 


Cm-1*—- 


: 1 [A««.W  d9  f ■>»] . 


we  can  get 


As  onr  original  values  of  Pr*.  and  P,%„  we  can  use  those  calculated 
for  determining  vortex  density  over  the  chord  fcy  the  law  of  the 


« 


jr(  J M3inci'\c<”r'.c03m8de'-*'.f  ?"Coj me'de^. 


DOC 


178  8 


PAGE  6 


cotangent  with  an  elliptical  distributing  function.  The  value  of  the 
angle  /,(«.*)  can  be  assuned  constant  (according  to  [3] 
according  to 


Proa  the  obtained  values  of  and  Awr<  ««  determine  according 

to  (1)  anfd  (6)  the  distributing  function  of  the  first  approxiaati  on 
and  also  P««  and  for  distribution  of  vortex  density  (6)  with  the 

variability  of  values  considered. 


He  know  that  the  most  intensive  formation  cf  free  vortices 
occurs  on  the  leading  lateral  edges.  The  angle  cf  slope  of  the 
vortices  to  the  lifting  surface  on  these  sections  can  be  determined 
approximately  in  the  fora  of 


l9Pj8''^  = 


) in  of 


« / 1 \ •(*)  f % 

coiof  ♦ (A  * yA  )P  CO)/ 

»«  \ •*»  ‘ >rv  «t«  l In 


17) 


n 


The  boundaries  of  the  sections  where  intensive  vortex  formation 
cccurs  are  very  difficult  to  define.  Moreover,  inside  the  boundary 
layer,  with  the  exception  of  the  wing  surface  itself,  the  velocities 
normal  to  it  are  generally  equal  to  zero  and  the  flow  lines  do  not 
adjoin  the  surface.  On  this  basis  we  assume  the  possibility  of 
directing  the  axes  of  the  free  vortices  along  these  flow  lines  and 
consider  the  angles  of  slope  of  the  free  vortices  not  equal  to  zero 


even  beyond  the  limits  of  the  section,  where  such  a branching  of  the 
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vortices  is  confirmed  e xper imental 1 y. 

If  we  assuie  that  the  distributing  functions  differs  little  from 
the  elliptical  and,  consequently,  the  deciding  role  in  the  formation 
of  the  downvash  is  played  by  vortices  of  great  intensity  flowing  off 
the  tip  sections,  then  there  can  be  no  great  error  in  "spreading" 
angle  8,  determined  by  expression  (7),  over  the  entire  wing  surface. 
This  assumption  does  not  contradict  the  condition  of  nonleakage, 
since  the  vortex  system  is  actually  formed  in  the  boundary  layer  and 
only  for  the  calculation  schematic  is  it  taken  to  the  middle  surface. 

The  second  and  subsequent  approximations  continue  until  the 
values  of  the  coefficients  of  the  series  coincide  within  the  limits 
of  assigned  calculation  accuracy. 


For  a small  number  of  breaks  in  the  leading  and  trailing  edges 
of  the  wing  the  calculation  can  be  reduced  to  solving  a system  of 
algebraic  equations,  which  has  definite  advantages  for  digital 
computer  calculation.  With  this  in  mind,  after  introducing  the 
symbols 


■jff 


H ’If)  . 


we  can  write  equation 
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J I I[r(n * ^ ^ . (8> 


where 


LM  L„. 


% . m 

r 

0 fu-t  ,n  0 /u*i,  n- 1 


If  we  substitute  in  (8)  the  value  ^ from  (6)  and  confine 
ourselves  to  the  number  of  q terms  of  the  series*  we  can  get  an 
equation  f or  (1  * 0 I,”*  of  coefficients  and  (m=  Given  q ♦ 1 

by  quantities  for  each  of  the  values  of  K*=  O”71,-1  (r we  can 
also  obtain  the  corresponding  number  of  algebraic  equations  which  are 
linear  with  respect  to  the  coefficients  of  series  (6), 


If  we  assuie  that  the  dovnwash  velocity  components  w et . are 
directed  in  planes  tangent  to  the  middle  surface  of  the  wing,  then 
for  the  component  of  total  relative  velocity  normal  to  the  axis  of 

the  vortex,  we  get 


V = 

It* 


=(co)a('+  Ww,)coji  - Wm3int  . 
y \ tt  trm  1 t«  c*«  r* 


In  this  case  the  coefficients  of  the  force  cf  the  pressure  of  the 
-th  section  normal  to  the  plane  of  the  projection  (direction  n,t,  in 
Pig,  2)  referred  to  an  element  of  the  projection  surface  Za„ d*',  is 
found  in  the  form  of 


(9) 
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however,  the  expressions  (c,)«.  and  (C*.)<«  are  found  by  multiplying  the 
integrand  function  (9)  by 


COJ  o( , 


jt  n 


and 


ctnfct  - d ) • respectively. 

V crtj 


Por  the  coefficient  cf  the  moment  of  hydrodynamic  forces 
relative  to  the  axis  running  through  the  leading  edge  of  the  root 
section  parallel  to  the  stand,  we  get 

If  we  assume  that  these  coefficients  are  constant  within  the  limits 
of  the  *,-th  band,  then  the  integral  characte r ist ic  can  be  obtained 

in  the  form  of 


^ 7 *«-•  1 

’ Z(v.*  KK  ' 


Us'.r  s,Yl, 


where  ar_,  represents  the  narrowing  at  the  beginning  of  the  £-th 
section,  - at  the  end. 


The  calculations  of  the  first  approximation,  performed  according 


- 


Figure  4 gives  as  an  idea  of  the  distributing  function  and  the 
distributed  hydrodynamic  characteristics. 
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Fig.  2.  Coordinate  system  of  wing  section.  KEY:  ( 1)  Middle  line  of 
^,-th  section,  (2)  Projection  plane. 
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Fig-  1-  Comparison  of  calculated  results  with  eiper inental . KEY:  (1) 
Calculation.  (2)  Sharp  edges,  (3)  Bound  and  elliptical  edges,  (4) 
Bartlett  and  Vidal,  (5)  Triangular  wing  X * 2- 
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